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abstract 

We study the asymptotic behavior of the B meson wave function in the framework of kx factorization 
- -. - , theorem. We first construct a definition of the /cT-dependent B meson wave function, which is free of 
' light-cone divergences. Next-to-leading-order corrections are then calculated based on this definition. The 
• treatment of different types of logarithms in the above corrections, including the Sudakov logarithms, and 
i those depending on a renormalization scale and on an infrared regulator, is summarized. The criticism 
Qh' raised in the literature on our resummation formalism and Sudakov effect is responded. We show that the B 
^ meson wave function remains normalizable after taking into account renormalization-group evolution effects, 



contrary to the observation derived in the collinear factorization theorem. 



o 



^ 1 INTRODUCTION 

O . 

[ The B meson distribution amplitude 0+(A;~^) plays an essential role in a perturbative analysis of exclusive B 
' meson decays based on collinear factorization theorem [1, 2, 3, 4, 5], where is the momentum carried by 
Cp ■ the light spectator quark. Its behavior certainly matters, and has been investigated in various approaches 
recently. Models of 0+ with an exponential tail in the large k^ region has been proposed [6]. Neglecting 
, ^ three-parton distribution amplitudes in a study by means of equations of motion [7, 8], was found to 

• be proportional to a step function with a sharp drop at large k~^ [9]. The asymptotic behavior of was 
^ , also extracted from an renormalization-group (RG) evolution equation derived in the framework of collinear 
Q I factorization theorem, which exhibits a decrease slower than l/k~^ [10]. That is, the B meson distribution 

• amplitude is not normalizable. This striking feature has been confirmed in a QCD sum rule analysis [11], 
J> ■ which includes next-to-leading-order (NLO) perturbative corrections. For a summary of the above progress, 

I refer to [12]. A similar divergent normalization of the B meson distribution function involved in inclusive 
■ decays has been observed recently [13]. 



A non-normalizable B meson distribution amplitude does not cause a problem in practice [14]. In a 
leading-order collinear factorization formula, only the first inverse moment X'^^{fi) = J dk~^(p^{k'^)/k^ is 
relevant [15, 16], which is a convergent quantity. The factor l/k~^ comes from a hard kernel of a decay mode. 
Note that a hard kernel would not be as simple as l/k~^ at higher orders, and information of more moments 
is required. However, the non-normalizability does introduce an ambiguity in defining the B meson decay 
constant Jb- The ambiguity can be understood through the matrix element, 

{0\q{y)Wy{n_)^Wo{n^)T h{{))\B{v)) = --^^{v ■ y, p^)tT (v /,„ , (1) 

where the coordinate of the anti-quark field q, y = (0, y~, Ot), is parallel to the null vector n_ = (0, 1, Or), 
h the rescaled b quark field characterized by the B meson velocity v, /i the renormalization scale, and T 
represents a Dirac matrix. The factor Wy(n„) denotes the Wilson line operator. 



Wy(n_) = Pexp 



roo 

ig / dXn_ ■ A{y + An_ 
Jo 



(2) 



The quantity -F(/u) is the HQET matrix element corresponding to the asymptotic value of the product 
Ibs/it^b in the heavy-quark limit. If the normalization (t)+{v • y = 0, /x) is divergent, the definition of fs 
demands a further arbitrary renormalization [17]. 
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We shall show that the above undesirable feature of 0+ is a consequence of adopting the collinear 
factorization theorem. It has been known that the collinear factorization formulas of many exclusive B 
meson decays suffer end-point singularities [18]. We regard these singularities as an indication [19] that 
the kx factorization theorem [20, 21, 22, 23, 24, 25] is more appropriate for studying these decays than 
the collinear factorization theorem. The perturbative QCD (PQCD) approach [26, 27, 28, 29] based on the 
kx factorization theorem has been developed. Retaining the parton transverse momenta kx [30], the end- 
point singularities disappear [31], and the resultant predictions are in agreement with most of experimental 
data [32]. Viewing these merits, it is very tempting to reanalyze the RG evolution effect on the B meson 
wave function (or the unintegrated B meson distribution amplitude) in the kx factorization theorem. Our 
conclusion is that the evolution effect does not drive the asymptotic behavior of the B meson wave function 
into 1/k^. Therefore, the B meson wave function is normalizable, and the B meson decay constant is 
well-defined. 

In Sec. II we find out a legitimate definition of the fcr-dependent B meson wave function. The NLO 

corrections to the B meson wave function are computed and compared to those in [10] in Sec. III. We 
respond to the criticism raised in [33, 34] on the PQCD formalism and on the Sudakov effect in Sect. IV. 
Section V is the conclusion. 

2 DEFINITIONS OF A WAVE FUNCTION 

We first construct the definition of the B meson wave function in the kx factorization theorem, which is 
nontrivial at all. Hence, our formalism for the NLO calculation differs from those in the literature [35, 36], 
which also involve the parton transverse degrees of freedom. A gauge-invariant definition of the B meson 
wave functions ^+{v ■ y, 6, /j.) is given via the nonlocal matrix element [25], 

{0\q{y)Wy{n.Yln_-y,oWo{n.) /i.Th{0)\B{v)) , (3) 

as a naive extension of Eq. (1) with y = (0,y~,b). The two Wilson lines Wy(n_) and Wo{n-) must be 
connected by a link In-;y,o at infinity in this case [25, 37]. 

As pointed out in [38], Eq. (3) contains additional collinear divergences from the region with a loop 
momentum parallel to n_. These light-cone divergences, cancelling each other as 6 — > 0, that is, as • 
y, h, /i) 4>+{v ■ y, fi), do not cause a problem in the collinear factorization theorem. In the k^ factorization 
theorem, however, they must be subtracted in a gauge-invariant way. Two modified definitions have been 
proposed in [38]: 

m{y)Wy{uyiu-y,oWo{u) jh.ThmB{v)) , (4) 

{0\q{y)Wy{n.)nn^;y,oWo{n_) jh.Th{0)\B{v)) 
{0\Wyin.)^Wyiu')In_-y,oll,y,oWoin-)Woiu'y\0) ' 

In Eq. (4) a non-light-like vector u has been substituted for the null vector n_, so that no collinear divergence 
is associated with the Wilson lines. In Eq. (5) n_ is maintained, but the light-cone divergences are regularized 
by the denominator, which contains the same light-cone divergences as in the numerator. As a gluon travels 
along n_, it does not resolve the detail of the valence quarks, which can then be replaced by the Wilson 
lines in an arbitrary direction u' [39, 40]. Both the above modifications with the appropriate Wilson links 
are gauge-invariant. Nevertheless, the universality of the B meson wave function is broken due to the 
appearance of the auxiliary scale, for example, = {k ■ u)/Vifi from Eq. (4). Fortunately, the evolution in 
the so-called Sudakov evolution [38], can be derived using the A;T-resummation technique [39, 41], such 
that the initial condition of the evolution remains universal. Note that Eqs. (4) and (5) do not approach the 
B meson distribution amplitude directly in the limit 6 ^ for general u and u', but convolutions of hard 
kernels with the B meson distribution amplitude [38, 42]. 

We have investigated the 0{as) diagrams in Fig. 1 according to both modifications, and found that 
Eq. (4) would change the ultraviolet structure of the quark- Wilson- line vertex correction in Eq. (3). This 
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problem can be explained using the pole term obtained from Fig. 1(c) (see Eq. (48) in the Appendix), 



Att 



T{e) 



mi 



(6) 



nig being an infrared regulator. If taking the u 



n_ 



I.e.. 



^ ^ oo limit in the above expression before 



making the expansion in e, only the first term contributes to the ultraviolet pole, which is 2/e in unit of 
C(sCf/{A7t), the same as in Eq. (3). If expanding the factor /mg)~'^ first (note that this expansion makes 
sense for ^ 0, i.e., C -h oo), the second term also contributes, and changes the ultraviolet pole into 1/e. 
Consequently, Eq. (4) is governed by a RG evolution different from that of Eq. (3). In this work we shall 
adopt Eq. (5), and demonstrate that the freedom in choosing the vector v! allows a correct RG evolution of 
the -B meson wave function. 



3 0{a,) CORRECTIONS 

The lowest-order evolution kernel for $+(A;"'", 6, /u) is written as 

ii:(°)(A;+, 6, n) = S{k+ - k'+) 



(7) 



which implies that the light spectator quark, carrying only a longitudinal momentum, is initially on-shell. 
It acquires the transverse degrees of freedom through collinear gluon exchanges, before participating a hard 
scattering [25]. As indicated in Eq. (7), we perform kx factorization in the conjugate b space. We then 
calculate the 0{as) corrections to Eq. (7) in dimensional regularization. A gluon mass nig is introduced 
to regularize the infrared divergences, so that we can clearly distinguish the ultraviolet poles 1/e from the 
infrared divergences represented by Inm^. As suggested in [38], a small plus component is added to the 
null vector ra_ in Eq. (5) at the intermediate step of the calculation. That is, we start with the Wilson line 
in a non-light-like direction u for the numerator, and take the u n_ limit eventually. Figures l(a)-l(g) 
contribute to the numerator of Eq. (5), and Figs. l(a)-l(d) with the quark lines being replaced by the Wilson 
lines along u' contribute to the denominator. 

To highlight the difference between the collinear and kx factorizations, we present the loop integrals 
associated with Figs. 1(a) and 1(b) and with Figs. 1(c) and 1(d), 



2e 



11 ■ V 



(27r)4-2e {vl + ie){P - mj + ie){u ■ I + ie) 
X [S{k+ - k'+) - (5 (fc+ - + 1+) exp(-zlT • b)] , 



2e 



(27r)4-2^ 



tr 



{k'-iy + ie 



75 /i- A+75 



m2 + ieu - 1 + ie 



[d{k+ - k'+) -5{k+ - k'+ + Z+) exp(-ilT • b) 



(8) 



(9) 



respectively, where the arbitrary Dirac matrix F has been set to 75. In the above integrals we have made 
explicit the ie prescription in the propagators l/v ■ I and 1/u ■ I, which follows the eikonal approximation of 
the quark or gluon propagators the loop momentum I flows through [23]. Note the additional Fourier factor 
exp(— ilT • b) associated with Figs. 1(b) and 1(d) [23, 41). In the collinear factorization theorem this Fourier 
factor disappears, corresponding to the 6^0 limit. Moreover, for u = n_, we obtain, from Eqs. (8) and 
(9), the counterterm identical to that found in the collinear factorization theorem [10], 



a.Cp 2 
Att e 



k+ e{k'+ - /,■+) e{k+ - k'+) 



+ 



k'+{k'+-k+)+ {k+-k'+) 



(10) 



The resultant anomalous dimension contributes to the splitting kernel in the RG evolution equation, that 
determines the asymptotic behavior of the B meson distribution amplitude [10]. The ultraviolet pole 1/e 
arises from the integration over the transverse loop momentum It up to infinity. It is then expected that 
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Eq. (10) will be absent in the kr factorization theorem due to the suppression in the large It region from 
exp(— ily • b). 

The 0{as) corrections N^^\k~^ , k'~^ ,b, /j,) from Figs. l(i), i = a---e, to the numerator of Eq. (5) are 
summarized below: 



.^ln(2.^) (l+ln^]5{k^-k'^), 



(11) 



OisCp 



\n{2u^)\n '\ S{k+-k'+) 



+4 



47r 

9ik'+ - fc+) 
{k'+ - k+)+ 



Ko {2u{k'+ - k+)b) - 4 



e{k+ - k'+) 
{k+ - k'+)+ 



Ait 



--21n2 + 21n + In ^ 



m„ 



m„ 



(12) 
(13) 



In 



In 



mo 



+ 



k+e{k'+ - k+) 
~ k'+{k'+ - k+)+ 

OisCp 



An 



-A 



k+e{k'+ ~ k+) 

k'+{k'+ - k+)+ 



6{k+ - k'+) 
Ko ( ^k+/k'+mgh ) - Kq {2u{k'+ - k+)b) 

S{k+ - k'+) 
Ko ( ^k+/k'+mgh \ - Kq {y2{k'+ - k+)h^ 



(14) 



In In —2 \ 

m„ A 3 



(15) 



For their detailed derivation, refer to the Appendix. The auxiliary parameter u = (n+ ■ u)/vv?, defined 
via C, = vk'^ , denotes the u dependence, is the Euler constant, and the subscript "+" in the factor 
l/{k'^ — k~^)+ represents the "plus" distribution. The self-energy corrections to the heavy quark field h and 
to the light spectator quark q are 



(1) _ 



at: 



9 




1 Airn^ ' 

- +ln^V 



S{k+ - k'+) , 



An 



1 + ilni^ 

2e 2 m^eTB 



- Sik+ - k'+) 



(16) 
(17) 



Some remarks are in order. The logarithms ln(2z/2) denote the light-cone collinear divergences mentioned 
before. Equation (11) does not contain a double pole l/e^ observed in [10] due to the replacement of the 
null vector n_ by the non-light-like vector u. For a similar reason, the single-pole term ln(^/fc+)/e, which 
leads to the type of Sudakov logarithms in the collinear factorization theorem [10], does not exist. As 
expected, Eq. (12), with the suppression from the Fourier factor, docs not generate the ultraviolet pole in 
Eq. (10). Because of the Bessel function Kq, the splitting effect from the plus distribution is negligible in 
the asymptotic region with large k'^. Equation (13) produces the double logarithm ln^(fc"'"/mg), which is 
not yet in the form of Sudakov logarithms in the fcy factorization theorem. After combining Eqs. (13) and 
(14), we derive the standard fcr-dependent infrared-finite Sudakov logarithms. 



27r 



ln2(A;+6)-(l-27E)ln(A;+6) 



(18) 



with the first (second) term being leading (next-to- leading). The u dependence is not made explicit, since 
it will be cancelled by that from the denominator. These logarithms should be resummed to all orders 
using the technique in [26]. The double infrared logarithm In^m^ in Eq. (15) is new, which does not 
exist in radiative corrections to a light meson process [35]. The important splitting effects, proportional to 
Ko{-sJk+ /k'+mgh) ^ \n{mgb), cancel between Eqs. (14) and (15). 
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The 0{as) corrections to the denominator are computed in a similar way, but with S {k'^ — k''^) being 
substituted for 5 (k^ — k'^ + /"'"). This substitution is made in that the denominator is to remove the hght- 
cone divergences arising from 0. Hence, the sphtting terms, i.e., the plus distributions in Eqs. (11)-(15), 
disappear. We choose u' = v for the incoming Wilson line (along the b quark), and a different u' for the 
outgoing Wilson line, such that the ultraviolet structure of the quark- Wilson-line vertex correction the same 
as in Eq. (3) is recovered. We emphasize that other choices of u' are equivalent, in view that the resultant 
B meson wave functions all collect the same soft structure of an exclusive decay. The expressions are 
summarized below: 



D 



(1) 



Ar(i) 
_asCF 

agCp 
An 



\u{2v^)\n 6{k+-k'+) , 



d{k+ - k'+) , 



\n{Au^iy 2) In ^— 6{k+ - , 



(19) 
(20) 

(21) 
(22) 



with the auxiliary parameter = (u' ■ n_)/vu'^. It is easy to check that the sum of the above corrections 
is free of the infrared cutoff mg. That is, the denominator in Eq. (5) does not alter the soft structure of the 
numerator, i.e., of Eq. (3), as requested above. According to our prescription, we set ln(4i/ ^) to unity. 
The total one-loop correction K^^^ to Eq. (7) is then written as. 



(1) 



j=a 



j=a 



OsCp 

Att 



+ In 



- + In (ire'^^ 11%'^^ - 2 ln'^{uk+b) + 2(1 - 2je) Hi^k+b) 



-(5 -27b) In 



nigb 



+ 2 In 



k+^b . rrinbe"^ 



In- 



+ 4-6^ 



37E 



(23) 



The ultraviolet pole 5/(2e) in unit of asCp/iA-K) is the same as the corresponding one derived in Eq. (8) 
of [10] under our prescription for fixing u' . Note that it differs from 3/e in [26], since it is the b quark 
field, instead of the rescaled one, that was adopted to define the B meson wave function in [26]. The pole 
5/(2e) should be partitioned in the way that 3/(2€) contributes to the factor F{ii) in Eq. (1) and 1/e to 
$+(A;+,6, /x). Here we do not perform such a partition. The splitting terms, which either vanish in the 
asymptotic region with large k'^ or cancel between Eqs. (14) and (15), have been dropped. 

The treatment of each term in Eq. (23) is explained as follows. The ultraviolet pole together with the 
constants are subtracted in a renormalization procedure. The logarithm ln(//6) is then summed to all orders 
using a standard RG evolution equation [10], giving an exponential R{b,fi,u). The Sudakov logarithms 
asln'^{iyk~^b) and ln(z^A;"'"6) are resummed, leading to the Sudakov factor S{k'^,b,u) [26]. The evolution 
of the B meson wave function from Eq. (5) is then given by 

^+ik+,b,^i) = Sik+,b,l^)Mk^A^i) , <^+(A;+,6,/x) = i?(6,/x,z.)0+(fe+,6,M = 1/6) ■ (24) 

The logarithms ln{mgb) and hi{k~^'^b/mg)ln{m,gb), representing the soft structure of the B meson wave 
function, are absorbed into the initial condition (p^{k'^ ,b, fi = 1/6) of the above evolution. They are then 
used to subtract the infrared divergences in the evaluation of hard kernels, i.e., in the so-called "matching" 
procedure. 

The exponentials in Eq. (24) are quoted from [26, 39] as 

'"'''^ djl 



S{k+,b, v) 
R{b,i^,u) 



exp 



exp 




In'— A{as{fi)) + B{u, a,(/i)) 



(25) 
(26) 
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with the one-loop anomalous dimensions, 



A 



(27) 



TT 



B = 



1 = 




^C7Ff5 + 21nz/2) 



(28) 



(29) 



where the running of the coupling constant as has been taken into account. It can be confirmed trivially that 
the exponent of Eq. (25) is identical to the Sudakov logarithms in Eq. (23), if the running of is frozen. In 
a practical analysis, the scale /x is set to a hard scale, which is usually of order The i/-dependences then 
cancel between S{k~^,b,i') and R(b,fj, = k~^,i'), such that the B meson wave function $+(fc^,6, /x = k~^) 
does not depend u. This cancellation is equivalent to that of the light-cone divergences, in agreement with 
the speculation in [38]. After the above cancellation, the Sudakov exponent in Eq. (25) reduces to the 
logarithms in Eq. (18), if neglecting the running of a^, and the anomalous dimension 7 is equal to —5 in 
unit of Q!s/(47r). We then bring the Wilson line direction u back to the null vector n_ as stated before. 

At last, we discuss the normalization of the B meson wave function $_|_(/c"'", 6, //) in the kx factorization 
theorem, which is defined as 



The Sudakov factor in Eq. (24) becomes identity in the limit b ^/k^, which approaches zero in the 
heavy quark limit for a fixed momentum fraction x = k~^ /{mBV~^). In the above limit the splitting terms 
proportional to the Bessel function Kq remain finite, and will not contribute to the evolution kernel. It is 
then obvious from Eq. (24) that the normalizability of the B meson distribution amplitude is not spoiled 
by the RG evolution effect, when evaluated according to Eq. (30). 



After completing the NLO calculation for the B meson wave function in the fcr factorization theorem, we are 
ready to respond to the cricitism raised by Dcscotes-Genon and Sachrajda [33] and by Lange and Neubert 
[34], which concerned the PQCD formalism and the Sudakov effect. We fully recognize that it is not easy 
to appreciate the delicacies of different approaches, and that misunderstandings are unavoidable. We hope 
that this section helps clarify these misunderstandings. 

It was concluded that a heavy-to-light form factor is not calculable in practice due to the ignorance of 
the heavy meson wave functions [33]. The word "calculable" is perhaps confusing. It is more appropriate to 
use "factorizable" , which means that a physical quantity can be written to all orders of as a factorization 
formula containing a hard kernel (Wilson coefficient) and wave functions, plus jet functions, Sudakov factors 
• • •. Then a form factor is factorizable in the PQCD approach based on the fc^ factorization theorem because 
of the absence of the end-point singularities. A form factor is not factorizable in QCD-improved factorization 
(QCDF) (considering only the leading contribution) [15], and partially factorizable in soft-collinear effective 
theory (SCET) [43], since both factorizable and nonfactorizable pieces exist at leading level. Therefore, a 
heavy meson wave function plays the role of an input in the PQCD approach the same as the form factor 
does in QCDF, and is determined by the value of the form factor from experimental data, lattice QCD, or 
sum rules. The inappropriate criticism in [33] is thus a conceptual misunderstanding of the PQCD approach. 

A difference between the pion form factor and the B ^ tt transition form factor was pointed out [33]: 
the former does not contain an end-point singularity in the collinear factorization theorem, but the latter 
does. Hence, it was questioned whether the PQCD approach, working for the former, can be extended to 
the latter. Our opinion is that both collinear and fcr factorization theorems are applicable to the pion form 
factor or to the decay B ^ jlv [44], and the numerical outcomes are not very different, because there is 
no end-point singularity. However, the end-point singularity in the collinear factorization formula of the 
B ^ TT form factor demands the use of the fc^ factorization theorem, which is more conservative than the 




(30) 



4 RESPONSE TO THE CRITICISM 
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collinear one: the parton transverse momenta should not be treated as a pure higher-power effect, when 

the end-point region of a parton momentum fraction is important. This is the motivation to develop the 
PQCD approach, and it makes sense to compare its predictions for two-body nonleptonic B meson decays 
with experimental data. 

It was stated that the b quark line a collinear gluon attaches could not be approximated by a Wilson line 
in the direction v, that is, could not be replaced by the rcscalcd b quark line [33]. In fact, the approximation 
holds, no matter a soft or collinear gluon attaches the b quark, in that it does not change the soft or collinear 
divergence of a loop integral. This is exactly the idea to remove the light-cone divergences in the numerator 
of Eq. (5) by introducing the denominator, where a Wilson line is substituted for the b quark. 

It has been argued that Fig. 1(c) does not contain the double logarithm lia^ (k ■ v / Vv^) , and the invariant 
{k ■ v)'^ /v'^ is irrelevant in the Sudakov resummation [26]. This argument was doubted in [33]. The explicit 
0{as) result in Eq. (15) has clarified the issue: ( is the only relevant invariant. It is well-known that a double 
logarithm arises from a vertex correction [45], such as Fig. 1(c), instead of from the type of corrections like 
Fig. 1(e). 

An expression equivalent to Eq. (25) was given by 
S{k'^, b, u) = exp 



( rC2i'k+ 




\ Jci/b 





In (|^^ j A{Ci, asifi)) + B{Ci,C2, a^ifl)) | , (31) 



where A{Ci,as) is equal to A{as) in Eq. (27) at one-loop level, and 

'g27B-lcf 



B{Ci,C2,as) = '^Cpln 



CI 



(32) 



Our choice Ci = C2 = 1 [26, 28], questioned in [33], is now justified, since it indeed reproduces the logarithms 
in Eq. (23). Other choices of Ci and C2 are equally fine: they lead to a change only in the next-to-leading 
logarithms, which can be compensated by the corresponding change in hard kernels. Hence, the choice of 
Ci and C2 is not a problem from the viewpoint of factorization theorem. 

It was claimed that the wave function in Eq. (4) defined in terms of the non-light-like vector u could 
not be related to the standard definition in terms of the null vector n_ [33]. As explained in Sec. II, the 
v dependence can be grouped into the Sudakov factor, and the initial condition of the Sudakov evolution 
is identified as the standard definition. The validity of Eq. (24) in the large b region was also challenged, 
because it suffers a large 0{as{l/b)) correction. The treatment of this correction, not multiplied by a 
logarithm, is a matter of factorization scheme: it corresponds to the constant terms in Eq. (23), and is 
allowed to shift freely between a wave function and a hard kernel. This shift is similar to that of the next-to- 
leading Sudakov logarithms resulted in by varying the parameters Ci^2 in Eq. (31). Therefore, it is always 
possible to choose a factorization scheme for a NLO evaluation of a hard kernel, such that Eq. (24) holds. 

The explicit ly dependence of the Sudakov factor derived from Eq. (4) [26] was pointed out in [33]. It 
has been known that this v dependence is cancelled by that of a soft function, which collects irreducible soft 
gluons to all orders. The cancellation has been demonstrated in the processes including Landshoff scattering 
[23], deep inelastic scattering [39], Drell-Yan production [39], inclusive semileptonic B meson decays [39], 
dijet production [46], and the B Dtt decays [47]. We believe that such a cancellation is general, though 
having not yet explored all other processes, since a physical quantity should not depend on this artificial 
dependence. 

There are two leading-twist B meson wave functions and [6], the latter being defined by, for 
example, the matrix element in Eq. (5) with /ri- replaced by /n^, where = (1,0, Oy) is another null 
vector. It was claimed [33] that the equality = $_ was assumed in the PQCD approach [28]. We make 
clear that this equality was never postulated. Precisely speaking, the B meson wave function adopted 
in [28] is identified as $_|- discussed here, and another wave function appearing as the combination 

— is numerically negligible as confirmed in [48]. The above combination was found to be important, 
when its contribution to a single term in the full expression of the form factor was investigated [33]. An 
observation based on such an incomplete analysis is certainly not solid. Because the wave function 
does not appear in the leading PQCD formalism, it is not urgent to discuss the corresponding Sudakov 
resummation. 
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It was concluded [34] that the Sudakov effect is not important for the soft contribution to the B ^ it 
form factor. First, we emphasize that the Sudakov logarithms studied in SCET [34] differ from what 
we discussed here and adopted in the PQCD approach [28]: the former appear in the Wilson coefficient 
associated with the soft form factor in the collinear factorization theorem, while the latter come from the 
wave functions in the kx factorization theorem. The difference manifests itself in the explicit expressions of 
the Sudakov factors: the latter is A;r-dcpcndcnt, but the former is not. Second, there is no conflict between 
the conclusions in [28] and in [34]. The weak Sudakov suppression in [34] refers to that on the whole form 
factor. The strong suppression in PQCD applies only to the end-point region of a momentum fraction (a 
form factor is factorizable in PQCD), and the suppression away from the end point is weak. Note that 
the strong Sudakov effect has been confirmed in [33] (see Page 271) for the model of the B meson wave 
function proposed in [28, 49]. Speaking of the whole form factor, whose contribution mainly arises from 
the non-end-point region, the suppression studied in PQCD is not significant either. For a more detailed 
explanation on this issue, refer to [50]. 

5 CONCLUSION 

In this paper we have surveyed the definitions of a wave function in the Ut factorization theorem given in 
Eqs. (3)- (5). The naive one in Eq. (3) contains additional light-cone collinear divergences, which cancel in a 
distribution amplitude in the collinear factorization theorem. These light-cone divergences are removed in 
the modified definitions of Eqs. (4) and (5) in a gauge-invariant way. However, Eq. (4), in which the Wilson 
line has been rotated away from the light cone to an arbitrary direction u, alters the ultraviolet structure 
of Eq. (3). Certainly, this change is not a problem, similar to the fact that the ultraviolet structure of a 
heavy-light currnt is changed under the HQET approximation. All the definitions of the B meson wave 
function are equivalent, as long as they collect the same soft structure of an exclusive decay. We have found 
that it is possible to maintain the ultraviolet structure by adopting Eq. (5). The dependence on a general 
u or u' can be factored into the Sudakov factor, such that the wave function, as the initial condition of the 
Sudakov evolution, is gauge-invariant and universal [38]. Eventually, the u or u' dependence of the Sudakov 
factor will be cancelled by that of a soft function as making predictions for a physical quantity. 

We have calculated the 0{as) corrections from Figs. 1(a)- 1(g) to the B meson wave function following 
the definition in Eq. (5), which contain three types of logarithms ln(A;"'"6), In(/Lt6) and \n{mgh). The leading 
and next-to-leading infrared-finite Sudakov logarithms ln(A;+6) have been verified, which are consistent with 
the Sudakov exponent adopted in our previous works. It has been observed that Figs. 1(b) and 1(d) do not 
generate the ultraviolet poles from the integration over the transverse loop momenta It due to the suppression 
from the Fourier factor exp(— zIt • b). Hence, the RG evolution from the summation of ln(/x6) is trivial. We 
have explained that the small h limit should be taken as 6 ^ 1 //c"*" in the factorization theorem, under 
which the Sudakov evolution factor becomes identity, and Figs. 1(b), 1(d) and 1(e) remain ultraviolet finite. 
As a consequence, the RG evolution effect does not spoil the normalizability of the B meson wave function, 
when evaluated according to Eq. (30). This is another indication that the fcr factorization theorem is a 
more appropriate framework for studying exclusive B meson decays than the collinear factorization theorem. 
At last, the infrared logarithms ln(mg6) are used to subtract the corresponding infrared divergences in the 
computation of hard kernels. 

Our NLO calculation is similar to that performed in [35], where the conjugate h space was also introduced. 
However, it was the 7*7 — > tt^ amplitude, instead of the pion wave function, that was studied in [35]. 
Therefore, the issues of the undesirable light-cone collinear divergences and of a legitimate definition of a 
A:T-dependent wave function were not addressed. The 0{as) corrections to the B — > ^Iv decay amplitude 
were computed in a different way in [36]. First, the issues mentioned above were not addressed either. 
Second, a parton was assumed to carry a transverse momentum initially, and the conjugate b space was not 
introduced. Third, a different type of double logarithms \v?{'mB/k^) was resummed, leading to the so-called 
threshold resummation [51]. Our opinion for proceeding a NLO analysis in the kx factorization theorem is 
that one must define a valid /cT-dependent wave function first (under a factorization scheme as stated in 
Sec. IV). Next, one computes the 0{as) corrections to the full parton-level amplitude, from which the wave 
function also evaluated at 0{as) is subtracted. This subtraction results in an infrared-finite hard kernel. 
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which is then substituted into a kr factorization formula to estimate the NLO effect. Therefore, our work 
provides a basis of the above systematic procedure. 
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National Science Council of R.O.C. under Grant No. NSC-92-2112-M-001-030 and by Taipei branch of the 
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A DETAIL OF CALCULATION 



We present the details of the 0{as) calculation in this Appendix. We assume the small plus component 
added to n_ to be negative, n"*" < 0, for convenience. One can always work out a loop integral, whose result 
is a function of u^, in the <0 {u^ > 0) region, and then analytically continue the result into the > 
{u^ < 0) region. Applying contour integration in the light-cone coordinates, we obtain the integral for the 
numerator of Eq. (5) associated with Fig. 1(a), 



u^v —uv^J (2TTr Jo 



u ' 



l+[2u-l+^ + u+{l^ + mj)] 



+ 



J—c 



l+[2v-l+'^ + v+{ll + ml)\ 



>6{k+ - k'+) 



(33) 



where the first term corresponds to the pole I = —u l'^ /u^ + ie for l'^ > 0, and the second term to the 
pole l~ = —v~l'^/v'^ — ie for Z+ < 0. The integration over It leads to 



OisCp 



U ■ V 



f 

Jo 



2tt u v 
°° dt 



mf, 



tl-25 



6{k+ - k'+) 



(34) 



where the variable change = nigt has been applied. It is easy to observe that the soft divergences in 
the above two terms cancel. Hence, the small parameter (5, introduced for convenience, will approach zero 
eventually. Working out the integration over i, we have 



agCp u ■ V I 47r/x 



T{e)B{6,€-6) 



2v- 



2u- 



d{k+ - k'+) , (35) 



whose 5 — > limit is given by 



u ■ V 



47r [u ■ v)"^ — v?v'^ 



""^'^ r(e)lnV^^(fc+-fc'+) 



(36) 



The above expression can be further simphfied into Eq. (11) as 0. 

Following the reasoning for Eq. (33), the loop integral associated with Fig. 1(b) is written as 



(1) 



2TTiu ■ V 



dHr 



r d^i 
J W 



u 



+e{k'+ - k+) exp(-ilT ■ b) 



U-V+ - u+v- J (27r)4 1 {k'+ - k+)[2u-{k'+ - k+f + u+{ll + ml)] 



v+0{k+ - k'+) exp(-zlT • b) | 
{k+ - k'+)[2v-{k+ - k'+Y + v+{ll + mf)] J ' 



(37) 
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where the ^-functions 9{k'^ — k'^) and 9{k'^ — k'^) correspond to the integration ranges Z+ > and Z+ < 0, 
respectively. The integration over It gives 



iV, 



(1) 



asCp u-v \ 9{k'+ -k+) ^^ 



l2u- 



ik'+ -k+f + mlb 



0{k+ - k'+) 
k+-k'+ 



'2v- 



ik+-k'+)^ + mp 



(38) 



We spht the above expression into 

Nl'\k+,k'+,b,i,) = 5{k+ - k'+) rdyN^'\k+,y,b,i,) + Nl'^{k+,k'+,b,iJ,) , (39) 

t/ 

where the first term can be combined with Na^\ and the second term defines the "plus" distribution. The 
first term is rewritten, by applying the variable changes y = {t + l)k'^ for the first Bessel function and 
y = (1 — t)k'^ for the second Bessel function, as 



rdyN^'\k+,y,b,fx) 
Jo 



agCp 



u ■ V 



-jlr^sKo{^2k^'t^ + mlb)Y 



(40) 



where the denominators t have been replaced by "^^ as in Eq. (34). In the asymptotic region with large 
A;"*", it is legitimate to extend the upper bound of t in the second integral from 1 to oo. Using the relation 



Eq. (40) becomes 



rdyNl'\k+,y,b,fi) = -^ 

Jo TT 



asCp u ■ V .,-1+5 



=2-'+'T{6)Ks{mgb) 



m„ 



f ^9 



4CbJ \2k+^b 



\J {u ■ vY — u^v'-^ 

Taking the limit (5 — 0, the 1/5 poles cancel between the above two integrals, and Eq. (42) reduces to 

u ■ V 



(41) 



(42) 



'dyNi'\k\yA^^) = -'"'^' 



2tt ^J(u ■ vY — u^v'^ w^v 



In — ■ In — ^- 



(43) 



The second term in Eq. (39) is written as 



u ■ V 



TT (u ■ vY — u'^V 

^{k+ - k'+) 



9{k'+ - k+ 
{k'+ - k+). 



-Ko {2u{k'+ - k+)b) 



{k+ - k'+)+ 



Ko (v^(A;+ - A;'+)6) 



(44) 



where the infrared regulators nig have been dropped, since a plus distribution is infrared finite. The com- 
bination of Eqs. (43) and (44) then gives Eq. (12). 

The loop integral associated with Fig. 1(c) from Eq. (9) is rewritten, in the light-cone coordinates, as 



2e 



2{k'+ - l+)u- 



(27r)4-2e [2(1+ - k'+)l- - ll]{2l+l- -ll- ml){u-l+ + u+t 



-(5(fe+ - k'+) . (45) 



10 



For < Z+ < we enclose the contour in the l~ plane over the pole l~ = {Ij, + mg)/(2Z+) — ie. For 

k'~^ < l^, we enclose the contour over the pole l~ = —u~l~^/u^ + ie. For I"*" < 0, there is no pinch singularity 
(noticing < in our choice), and the contour integration vanishes. Hence, Eq. (45) becomes 



dr 



{k'+ - l+)l+u- 



ik 



dl+- 



(27r)4-2e I Jo W+ll + {k'+ - l+)m1][u+{ll + m^) + 2u-l+^] 

(/+ - k'+)u+u- ' 



k'+ \u 



+ 2u-l+{l+ - k'+)][u+{ll + ml) + 2«-Z+2] 



\6{k+ -k'+) . 



(46) 



The integration over It leads to 

agCp ( 47r/Lt^ 



27r I 



rfe) 



dtt 



-1+25/ 



+ 



'AC 



mi 



dtt 



-2e 



dtt-^{t -ly 



dtt-^+^' ( Kt'' + 1' 



mi 



(47) 



where the variable change = k'~^t has been made. The above expression has been arranged in a way that 
the infrared divergences from t ^ cancel in the first line, and the linear ultraviolet divergences cancel in 
the second line. We have 



2tt I m^ 



Fie) { 3(26,1 -e) 



2\mj 



+ 



'AC 



B{1 - e, 2e) + ( ^ ) 5(1 - e, 2e) - ^ \5{k+ - k'+) , (48) 



mi 



l-2e 

\ a / 

which leads to Eq. (13) by employing the expansion, 

1 



r(e) 



1 - 7Ee + 



2 12 i 



(49) 



We calculate the correction from Fig. 1(d) in the light-cone coordinates, 

dH 2{k'+ - l+)u- 



N. 



(1) 



(27r)4 [2{l+ - k'+)l^ - ll]{2in- -I'^T- m^g){u-l+ + u+t 
x6{k+ - k'+ + /+) exp(-ilT • b) . 



(50) 



Applying the reasoning for Eq. (45) , we have 



N 



(1) 



-ig^CFA-Kieik'^ - k+) J ^ 



k+{k'+ - k+)u- exp(-ar • b) 



(27r)4 \ [k'+l^ + k+mj][u+{ll + mf, + 2u-{k'+ - k+f 
k'^u^u~ exp(— zlr ■ b) 1 



[u+ll - 2u-k+{k'+ - k+)][u+{ll + ml) + 2u-{k'+ - k+)^] 
The integration over It gives 

(1) _ c(sCf k+9{k'+ - k+) / / .\ „ / I2u- 



(51) 



N. 



TT k'+{k'+ - k+) 



Ko Jk+/k'+mgb - Ko 



u 



+ 



{k'+ - k+)^ + m% 



(52) 
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We then adopt the sphtting similar to Eq. (39), whose first term is rewritten, by applying the variable 
change y = — t) for the first Bessel function and y = {t + for the second Bessel function, as 



dyN^^\k+,y, b, ii) = dtt-^+^^Ko {VT^nigb) - dtt-^+^^Ko (^^H^ + 



oo 

+ 1 -^Ko{2Cbt) 



/o t + 1 

It is easy to show, using the relation, 



—ax = — CSC (z/tt) 

x+1 2 ^ ^ 



ta) 



(53) 



(54) 



where Jo denotes the Anger function, that the third term is given by 7r/(4^6) in the asymptotic region with 
large k^. Employing Eq. (41), Eq. (53) becomes 



rdyN^'\k+,y,b,i,) 

J 



27r 



In 



In — ^ 1 

m„ 2 6 



(55) 



The second term in the splitting with the plus distribution can be obtained in a way similar to Eq. (44). 
The 0{a.s) correction from Fig. 1(e) is written as 



i 2^ f dH 



n^'- /) 



{k'-l)'^{P-mj) 
x6{k+- k'+ + /+) exp(-zlr • b) , 



75 jh- /i+75 



vl 



dH 



2{k'+ - l+)v- 



(27r)4 [2{l+ - k'+)l~ - l^]{2l+l- -ll- ml){v-l+ + v+l-) 



'.5{k+ - k'+ + 1+) exp(-ilT • b) . 



(56) 



The substitution of v for n introduces an additional pole = v^l^ /v^ — ie in the region < k'~^ . It is 
straightforward to confirm that the contribution from this additional pole vanishes. Hence, the result of 
Fig. 1(e) is similar to that of Fig. 1(d), but with the variable replaced by k ■ v/Vv^, which is Eq. (15). 
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(a) (6) (c) id) 




(e) if) (g) 



FIG. 1 0{as) diagrams for the B meson wave function. 
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